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(57) ABSTRACT 

A system and method of measuring a residence time in a 
gas-turbine engine is provided, whereby the method includes 
placing pressure sensors at a combustor entrance and at a 
turbine exit of the gas-turbine engine and measuring a com- 
bustor pressure at the combustor entrance and a turbine exit 
pressure at the turbine exit. The method further includes 
computing cross-spectrum functions between a combustor 
pressure sensor signal from the measured combustorpressure 
and a turbine exit pressure sensor signal from the measured 
turbine exit pressure, applying a linear curve fit to the cross- 
spectrum functions, and computing a post-combustion resi- 
dence time from the linear curve fit. 
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METHOD FOR MAKING MEASUREMENTS 
OF THE POST-COMBUSTION RESIDENCE 
TIME IN A GAS TURBINE ENGINE 

CROSS-REFERENCE TO RELATED 5 

APPLICATIONS 

This application claims the benefit of U.S. Provisional 
Patent application Ser. No. 61/503,855 entitled “METHOD 
FOR MAKING MEASUREMENTS OF THE POST-COM- 
BUSTION RESIDENCE TIME IN A GAS TURBINE 
ENGINE” filed on Jul. 1, 2011. The entirety of the above- 
noted application is incorporated by reference herein. 

ORIGIN OF THE INVENTION 15 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used only by or for the Government for Government 
purposes without the payment of any royalties thereon or 20 
therefore. 

BACKGROUND 

A challenging issue confronting the air transport system is 25 
the demand for the reduction of the emissions of oxides of 
nitrogen. The formation of thermal NO,^ in a gas turbine 
engine depends on the stoichiometry, the residence time lin- 
early, and on the reaction temperature exponentially. Zeldov- 
ich thermal NO,^ may be produced by oxidation of atmo- 30 
spheric nitrogen in post flame gases. As turbine blade 
resistance to high temperatures improves, nitrogen produc- 
tion in the post-combustion zone may become more impor- 
tant. While residence time is not as signiflcant as temperature 
in formula predicting NO,^ production, it is a necessary factor 35 
and should be as accurate as possible. The characteristic 
combustor residence time can be defined as the ratio of the 
combustor volume to the bulk (volumetric) flow rate. This 
value is estimated from geometry and operational data. 
Detailed geometrical and operational data from gas turbine 40 
engine manufacturers, however, is frequently unavailable. 
Furthermore, post-combustion residence time measurements 
are not available to verify analytical estimates. Consequently, 
in order to improve the technology to satisfy future emission 
prediction goals, a different concept for determining the char- 45 
acteristic post-combustor residence time is required. 

SUMMARY 

The following presents a simplified summary in order to 50 
provide a basic understanding of some aspects of the innova- 
tion. This summary is not an extensive overview of the inno- 
vation. It is not intended to identify key/critical elements or to 
delineate the scope of the innovation. Its sole purpose is to 
present some concepts of the innovation in a simplified form 55 
as a prelude to the more detailed description that is presented 
later. 

In an aspect of the innovation a different concept for deter- 
mining characteristics of post-combustor residence time is 
disclosed. The concept is based on determining the post- 60 
combustion residence time delay due to convection of entropy 
at the flow speed in the combustor of a gas turbine engine over 
a range of operating conditions. It is evaluated from the slope 
of the combustor sensor signal/turbine exit signal pressure 
cross-spectrum phase angle over an appropriate frequency 65 
range where the measured signal can be attributed to indirect 
combustion noise. 
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In another aspect of the innovation the innovation, a 
method of measuring a residence time in a gas-turbine engine 
is provided, whereby the method includes placing pressure 
sensors at a combustor entrance and at a turbine exit of the 
gas-turbine engine and measuring a combustor pressure at the 
combustor entrance and a turbine exit pressure at the turbine 
exit. The method further includes computing cross-spectrum 
functions between a combustor pressure sensor signal from 
the measured combustor pressure and a turbine exit pressure 
sensor signal from the measured turbine exit pressure, apply- 
ing a linear curve fit to the cross-spectrum functions, and 
computing a post-combustion residence time from the linear 
curve fit. 

In yet another aspect of the innovation the innovation, a 
system to measure a post-combustion residence time in a 
gas-turbine engine is provided and includes a measurement 
component that measures a plurality of combustion and/or 
turbine flow pressures in the engine, a receiving/calculation 
component that calculates a plurality of cross-spectrum func- 
tions based on the measured plurality of combustion and/or 
turbine flow pressures, a tabulation component that tabulates 
a linear curve fit based on the calculated plurality of the 
cross-spectrum functions, and a computation component that 
computes the post-combustion residence time in the gas- 
mrbine engine. 

To accomplish the foregoing and related ends, certain illus- 
trative aspects of the innovation are described herein in con- 
nection with the following description and the annexed draw- 
ings. These aspects are indicative, however, of but a few of the 
various ways in wliich the principles of the innovation can be 
employed and the subject innovation is intended to include all 
such aspects and their equivalents. Other advantages and 
novel features of the innovation will become apparent from 
the following detailed description of the innovation when 
considered in conjunction with the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. I is a perspective view of an example gas-turbine 
engine with instrumentation that can be utilized to measure 
post-combustion residence time in accordance with an aspect 
of the innovation. 

FIG. 2 is a schematic illustration of a turbine-combustor- 
tailpipe noise system diagram in accordance with an aspect of 
the innovation. 

FIGS. 3-9 are graphical illustrations of a magnitude 
squared aligned coherence (MSC) function in part (a) and a 
cross-spectrum phase angle in part (b) in accordance with an 
aspect of the innovation. 

FIGS. 10-16 are graphical illustrations of a magnitude 
squared aligned coherence (MSC) function in part (a) and a 
cross-spectrum phase angle in part (b) in accordance with an 
aspect of the innovation. 

FIG. 17 is a graphical illustration of post-combustion resi- 
dence times as a function of engine power in accordance with 
an aspect of the innovation. 

FIG. 18 illustrates an example system incorporating a 
method of measuring post-combustion residence time in a 
gas-turbine engine in accordance with an aspect of the inno- 
vation. 

FIG. 19 illustrates an example flow chart of a procedure 
that measures post-combustion residence time in a gas-tur- 
bine engine in accordance with an aspect of the innovation. 

FIG. 20 illustrates a block diagram of a computer operable 
to execute the disclosed architecture. 

DETAILED DESCRIPTION 

The iimovation is now described with reference to the 
drawings, wherein like reference numerals are used to refer to 
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like elements throughout. In the following description, for 
purposes of explanation, numerous specific details are set 
forth in order to provide a thorough understanding of the 
subject innovation. It may be evident, however, that the inno- 
vation can be practiced without these specific details. In other 5 
instances, well-known structures and devices are shown in 
block diagram form in order to facilitate describing the inno- 
vation. 

While specific characteristics are described herein (e.g., 
thickness), it is to be understood that the features, functions to 
and benefits of the innovation can employ characteristics that 
vary from those described herein. These alternatives are to be 
included within the scope of the innovation and claims 
appended hereto. 

While, for purposes of simplicity of explanation, the one or 15 
more methodologies shown herein, e.g., in the form of a flow 
chart, are shown and described as a series of acts, it is to be 
understood and appreciated that the subject innovation is not 
limited by the order of acts, as some acts may, in accordance 
with the innovation, occur in a different order and/or concur- 20 
rently with other acts from that shown and described herein. 

For example, those skilled in the art will understand and 
appreciate that a methodology could alternatively be repre- 
sented as a series of interrelated states or events, such as in a 
state diagram. Moreover, not all illustrated acts may be 25 
required to implement a methodology in accordance with the 
innovation. 

Reducing NO,^ emissions and aircraft fuel bum by more 
than 75% while achieving perceived cumulative noise levels 
7 1 dB below stage 4 limits in subsonic vehicles is a ftimre goal 30 
of the aerospace industry. Over the years studies have been 
conducted on gas turbine (turbofan) engines to evaluate mr- 
bine tone generation, attenuation of direct combustion noise, 
etc. Indirect combustion noise, on the other hand, was ini- 
tially thought to be a non-contributor to the engine core noise 35 
and, thus, was investigated only analytically and in model 
scale tests. 

Low frequency noise generated in the turbo fan engine core 
may make a significant contribution to the overall noise sig- 
nature in the aft direction at the low power settings, which are 40 
used on an airport flight approach trajectory. This type of low 
frequency noise may be a problem for future aircraft. Two 
possible low frequency noise sources are “direct” and “indi- 
rect” combustion noise. The source of combustion noise 
attributed to the unsteady pressures produced by the unsteady 45 
combustion process that propagate through the turbine to the 
far field is called the “direct” combustion noise source. The 
other source of turbofan engine combustion noise is known as 
the “indirect” mechanism in which the noise is generated in 
the turbine by the interaction of entropy fluctuations, which 50 
also originate from the unsteady combustion process, as they 
propagate through regions characterized by mean flow veloc- 
ity or pressure gradients in the turbine stages. The innovation 
disclosed herein uses measured data from the indirect com- 
bustion noise in the combustor and the turbine exit to directly 55 
measure the post-combustion residence time. 

The formation of thermal NO,^ depends on the stoichiom- 
etry, the residence time linearly, and on the reaction tempera- 
ture exponentially. Zeldovich thermal NO,^ may be produced 
by oxidation of atmospheric nitrogen in the post flame gases. 60 
As turbine blade resistance to high temperatures improves, 
nitrogen production in the post-combustion zone may 
become more important. NO„. levels increase with increasing 
post-combustion residence time. The characteristic combus- 
tor residence time can be defined as the ratio of the combustor 65 
volume to the bulk (volumetric) flow rate. This information, 
however, is proprietary, which limits the information avail- 
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able about the post combustion residence time in current and 
future combustor design concepts. Thus, the innovation dis- 
closes an alternate concept for determining the characteristic 
post-combustor residence time. The innovation determines 
the post-combustion residence time delay due to convection 
of entropy at the flow speed in the combustor. It is evaluated 
from the slope of the combustor sensor signal/turbine exit 
signal pressure cross-spectrum phase angle over an appropri- 
ate frequency range where the measured signal can be attrib- 
uted to indirect combustion noise. 

The net travel time of the indirect combustion noise signal 
from the combustor to the turbine exit and the far field 
increases since the travel velocity of the entropy fluctuations 
to the turbine is the flow velocity, which is a small fraction of 
the speed of sound, in the combustor. The innovation demon- 
strates that the pressure and entropy should be in phase in the 
combustor. As a result, the pressure signal from an indirect 
combustion noise source would be delayed relative to a pres- 
sure signal from a direct combustion noise source since an 
indirect combustion noise signal does not travel with the 
speed of an acoustic wave until it interacts with the turbine. 

The iimovation shows that the cross-spectra and correla- 
tion function between a combustor sensor and far-field micro- 
phones are tools that provide a way to separate “direct” and 
“indirect” coherent combustion noise due to this travel delay 
time. This paper uses measurements in the combustor and 
mrbine exit from a test engine to directly measure the post- 
combustion residence time. The innovation uses tools that are 
part of signal -processing theory to study a combustor pres- 
sure sensor signal and a turbine exit pressure sensor signal. 
The cross-spectral density phase measurement identifies a 
time delay that corresponds to the convective time delay. The 
magnitude of the coherence between the two sensors identi- 
fies the spectral region of importance as being in the 50-250 
Hz frequency range. The innovation discloses the cross-spec- 
tral density phase angle and the coherence over a range of 
operating conditions and calculates the post-combustion resi- 
dence time from the slope of the cross-spectral density phase 
angle. 

Prior to disclosing the innovation, the information dis- 
closed herein will be presented as follow: 1 ) First expressions 
for NO„. emission; 2) First expressions for residence time; 3) 
Engine noise data; 4) A linear system theory; 5) A system 
model; 6) Results; and 7) The post-combustion residence 
time results calculated from the cross-spectrum between a 
combustor pressure sensor and a turbine exit pressure sensor 
are presented. 

The operation of a gas turbine engine can be correlated 
with NO,^ emission levels using regression analysis of mea- 
surable test parameters or by consideration of time scales and 
chemical kinetics or using both sets of variables along with 
such variables as equivalence ratio, fuel flow rate and Mach 
number. A correlation of NO,^ emission levels determined by 
others for propane air combustion is based only on the resi- 
dence time and adiabatic flame temperature. It was deter- 
mined that the NO levels are principally a function of adia- 
batic flame temperature and combustor residence time, which 
is represented by the expression: 


Eno, = f„.exp(-72.28 + 2.80V 

where is the NO„. emission index (gN02/kg-fuel), 

'^adiabatic 1^6 adiabatic flame temperature (° K), and is 
the combustor residence time (ms). They found that over the 
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range of pressures from 5 to 30 atmospheres, there is no 
significant observed departure from this expression for inlet 
temperatures 727K and higher. 

Combustor residence time is given by the bulk 

flow through the combustor volume expressed in equation 5 

( 2 ): 


volume VP (2) 

volumetricflowrate mRT 10 


Others have stated that that the residence time in a conven- 
tional combustor and in a micro-combustor is approximately 
7 ms and 0.5 ms respectively. 

The primary combustion zone residence time can also be 
calculated by equation (3): 


^PZPcombair (3) 

TP7 = 2U 

>hpz Tf,D,PZ 

In examining a set of exemplified gas turbine dual fuel, dry 
low emission combustion system, primary zone residence 
times were found to be 2.71, 1.35, 8.17, 4.09, 9.84, and 4.92 25 


Gas turbine NO,^ production, however, is more complicated 
than NO,j emission from a propane combustor, described 
above, since in addition to the resident time dependence the 
reaction rate is assumed to be a function of pressure in addi- 30 
tion to temperature or: 

reaction rate=p"‘exp(z7) (4) 

and the mixing rates are assumed to be a function of linear 
pressure drop or: 35 

mixing rate=(AP/P)A; (5) 

Consequently, 


/ PV \/AP\^_ 

_ AVc[APfF^‘*"'~’‘^txp{zT) 
mT 


( 6 ) 


45 


Others have correlated a large set of engine data using A=9x 
10“^, x=0, m=0.25, and z=0.01 so that: 


^NOx 


9 .0 X 1 0“ V, pp^ exp(0.0 1 Tst ) 

^a'^PZ 


(7) 50 


Still others recast equation (7), thereby making changes to 
improve correlation with data and derived: 55 


doing a least squared curve fit to a large data set. The primary 
zone residence time formulation has become more complex 
as combustor design has become more complex. In part, this 
may be due to it not being a measurable quantity. When used 
as correlation factor it should be referred to as a primary zone 
NO,^ emission residence time and not as the primary zone 
residence time. As mentioned above, the innovation disclosed 
herein uses a procedure to measure the post-combustion resi- 
dence time using signal processing methods. As a conse- 
quence, tliis post-combustion/post-flame residence time 
becomes available for consideration with knowledge of any 
engine company proprietary combustor geometry design 
information or proprietary operating parameters. 

Referring now to the figures, FIG. 1 is an example engine 
100 that can be utilized to conduct experiments to obtain 
engine test data and FIG. 2 illustrates a turbine-combustor- 
tailpipe noise system diagram 200 in accordance with an 
aspect of the innovation. The engine 100 is a dual-spool, 
mrbofan engine that has a direct drive, a wide chord fan 
connected by a long shaft to a low-pressure turbine spool, and 
a high-pressure compressor connected by a concentric short 
shaft to a turbine high-pressure spool. The fan diameter is 
approximately 34.2 in. The combustor design is a straight- 
tlirough-flow annular geometry with 1 6 fuel nozzles and 2 
igniters. 

The engine-internal instrumentation in this configuration 
includes a higli-temperature pressure sensor with air cooling 
in a combustor igniter port 102 (hereinafter “combustor sen- 
sor”), a first 104 liigh-temperature pressure sensor with air 
cooling at a first turbine exit (hereinafter “first turbine exit 
sensor”) and a second 106 liigh-temperature pressure sensor 
with air cooling at a second turbine exit (hereinafter “second 
turbine exit sensor”). 

The data acquisition system had a sampling rate of 65 536 
Hz and a sampling duration of approximately 70 s. The spec- 
tra were calculated using a 50 percent overlap, which permit- 
ted data reduction using approximately 254 overlapped 
ensemble averages at a bandwidth resolution of 2 Hz. Signal 
estimation parameters are shown in Table 1 below. 


TABLE 1 


Soectral estimate parameters. 

Parameter 

Value 

Segment length, (data points per segment), N 

32 768 

Sample rate, r^, samples/s 

65 536 

Segment length, = N/r^, s 

0.500 

Sampling interval. At = l/r^, s 

1/65 536 

Bandwidth resolution, = Af = 1/T^ = r^/N, Hz 

2.0 

Upper frequency limit, = l/2At = tJ2, Hz 

32 768 

Propagation time delay/lag (T = 9° C., r = 

0.09044 

30.43 m) Tp = 5 927/65 536, s 


Number of independent samples, n^ 

128 

Overlap 

0.50 

Sample length, g 

«70 


1.5xl0‘hrjvo„ -0.5T„)‘’-^exp(-71100/7’s,) (8) 

~ F^oOiAPi/Pif-^ 

where is the residence time in the NOx production 
region, x^^ is the evaporation time, T^, is the reaction tempera- 
ture, ?3 is the combustor inlet pressure, and (AP3/P3) is the 
combustor pressure drop. 

Consequently, the primary zone residence time has evolved 
to become a NO„. emission production parameter evaluated by 


The spectra and cross-spectra are estimated using a non- 
parametric method, which is based on averaging multiple 
windowed periodograms using overlapping time sequences. 
60 Using these spectra and cross spectra, the magnitude squared 
coherence is calculated to measure the similarity of the ampli- 
tude variations at particular frequencies. The " accent will be 
used to denote the statistical basis of a variable. This is done 
to avoid confusion with calculations of coherence using a 
65 single segment or block which yield a coherence of unity. The 
concept used is based on determining the post-combustion 
residence time delay due to convection of entropy at the flow 
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speed in the combustor. It is evaluated from the slope of the 
combustor sensor signal/turbine exit signal pressure cross- 
spectrum phase angle over an appropriate frequency range 
where the measured signal can be attributed to indirect com- 
bustion noise. 

The appropriate frequency range is determined from the 
combustor sensor signal/turbine exit signal magnitude 
squared aligned coherence (MSC) function illustrated in 
equation (9) below. 


8 

The methods used herein are based on system theory devel- 
oped for linear systems with random inputs. The linear sys- 
tem theory disclosed herein is in the frequency domain. The 
output spectral density function, G and the cross-spectra den- 
5 sity function, G, is related to an input spectral density func- 
tion, G through frequency response function, H, representing 
the turbine as 

(11) 

m and 


^y,y 

FIGS. 3-9 illustrate the MSC function in part (a) and the 
cross-spectrum phase angle in part (b) where the post-com- 
bustion residence time is at various percentages of maximum 
power as measured by the combustor sensor 102 and the first 
turbine exit sensor 104. Specifically, FIGS. 3-9 are at 48, 54, 
60, 71, 87, 98, and 99 percent of maximum power respec- 
tively. 

Further, FIGS. 10-16 illustrate the MSC function in part (a) 
and the cross-spectrum phase angle in part (b) where the 
post-combustion residence time is at various percentages of 
maximum power as measured by the combustor sensor 102 
and the second turbine exit sensor 106. Specifically, FIGS. 
10-16 are at 48, 54, 60, 7 1 , 87, 98, and 99 percent of maximum 
power respectively. 

Also illustrated in part (a) of FIGS. 3-16 is a coherence 
threshold calculated from: 

(10) 

where this formula determines a P-percent threshold confi- 
dence interval using the number of data segments/blocks, n,,, 
used in the periodogram method spectral estimator. The 95 
percent confidence interval based on n=128 independent 
samples is 0.0233. The spectra are calculated using a 50 
percent overlap and the 95 percent confidence interval based 
onn=273 samples is 0.0109. These indicators show the MSC 
function is reliable up to about 400 Hz. However, MSC func- 
tion is above 0.1 only in a region from 30-250 Hz. The phase 
angle variation in this region is attributed to indirect combus- 
tion noise. 

Instead of relying on the confidence interval given by Eq. 
(10), which is based on a statistical theory, to obtain a thresh- 
old value for y„„^(n^), a deliberately unaligned time history 
can be used to create the threshold value. If one of the time 
histories is shifted by a time delay more than the segment/ 
block length, Td=N/r^i, then the two time histories are totally 
independent unless tones are present. This deliberate de-cor- 
relation establishes a coherence threshold and also identifies 
any tones in the signals. Shifting the signals by this time delay 
removes the coherence of random noise but leaves the coher- 
ence of periodic functions which are sometimes identified as 
hidden periodicities, concealed spectral lines, or un-damped 
sinusoids in noise. The deliberately unaligned coherence is 
also shown in part (a) of FIGS. 3-9. Note that the higher 
statistical confidence interval based on the number of inde- 
pendent records (n^=128) is a more conservative estimate of 
the measured coherence threshold. The statistical coherence 
threshold can be used with confidence since it can be com- 
pared with a measured coherence threshold. The coherence 
value is below the 95 percent statistical confidence interval 
above 400 Hz. Consequently, this is the upper frequency limit 
for which data is available for analysis using a linear system 
model. 


Gx^=Hxy<J)Gxxr ( 12 ) 

where x is the input signal from the high-temperature pres- 
sure sensor with air cooling in a combustor igniter port 102 
and y is the output signal from either the first 104 or second 
106 high-temperature pressure sensor with air cooling at the 
mrbine exit. 

The cross spectral density and the frequency response 
functions are complex valued quantities, which can be 
expressed in terms of a magnitude and an associated phase 
angle. Tliis will be expressed herein using complex polar 
notation. 


Gxyl/)=IGxyiy)leXp[/m] 

(13) 

Hxy(/)=IHxy(/)-exp[-/H>(/)] 

(14) 


Before plotting the cross spectral density phase angle, 
phase unwrapping is applied to the phase angle to avoid a 
jump of 2jt in the phase caused by the ATAN2 function. 

The system under consideration has a combustion noise 
input with a measured spectrum, G, which includes acoustic 
and hydrodynamic components. The system measured output 
quantities are assumed to be related as follows: 


W=-?^9,.o'”(fi*S9.9'"(/) 

(15) 

'”(fi=l^9,.o'"(/)l"*S9,9”(/) 

(16) 


and the measured MSC by: 


.-,.2 1 ^ 9,101 

Z9.10 “ xm xm 

*^ 9 , 9 *^ 10,10 

where m indicates noise may be included in the measured 
quantities. The unknown that will be identified is the turbine 
frequency response function, H"'(f) at a range of operating 
conditions. 

The system model disclosed herein is applied in the 50-250 
Hz frequency range. The model involves the turbine attenu- 
ation and the convective time delay of the of the entropy 
signal. The plant being modeled is the turbine. The input to 
the plant is the total pressure signal and the measurement 
made in the combustor is of the total pressure signal. Conse- 
quently, the available input auto-spectrum is that of the total 
pressure signal. To aid in physical interpretation, a standard 
template parametric model form will be used. The model 
form disclosed herein is in a parametric reduced order fre- 
quency domain representation. The parameters depend non- 
linearly on the operating point. However, at each operating 
condition the system will be assumed to be linear and the 
same parametric form will be used so that source separation 
will be obvious. The nonlinear operation is then described by 
a linear model at a range of observed operating points each 
identified by a set of parameters. As mentioned above, the 
iimovation discloses measurements at the 48, 54, 60, 87, 98, 
and 99 percent maximum power settings. Over the frequency 
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range 50-250 Hz, the turbine exit signal will be the result of an 
attenuation of the input signal by an amount K. At the turbine 
exit, the input signal will also have a time delay since the 
indirect combustion noise travels at the flow velocity t. 

Consider an input signal x(t) with a spectrum G,„(f) for a 
system with transfer function H,^^(f) and output signal y(t). 
Then the cross spectrum is given by Equation (12). For the 
turbofan engine 100, the input to the turbine is the direct 
acoustic signal x^(t) and the time delayed entropy signal, x,.(t), 
with a delay of The entropy signal may represent in addi- 
tion to a temperature fluctuation moving with the flow any 
other disturbance moving with the flow such as a vorticity 
fluctuation. 

Referring back to FIG. 2, as mentioned above, FIG. 2 
illustrates a turbine-combustor-tailpipe noise system diagram 
200 in accordance with an aspect of the innovation. The 
output turbine noise signals are y„/(t) and y,(t). The output 
signal, y, is the sum of the direct combustion noise signal, y^, 
and the indirect combustion noise signal, y,: 

y=yj^yt (i8) 

The direct combustion noise cross-spectral density, 
is a product of the direct combustion noise turbine transfer 
function, H^(f), and the direct combustion noise input spec- 
tral density, G,^^. 25 


^xdyd~^JJ)'^^xd 

The indirect combustion 


•^xdyi' 


15 


20 


(19) 

cross-spectral density. 


is a product of the indirect combustion noise turbine 
transfer function, H,(f), the time delay factor, and the 

input indirect combustion noise spectral density, G,^,., which 
corresponds to an equivalent fluctuating entropy spectral den- 
sity function. 

( 20 ) 

The indirect combustion noise turbine transfer fimction, 
H,(f), is assumed to have a representation, H(f). The direct 
combustion noise turbine transfer function, H^(f), is assumed 
to have a representation, a H(f). Where a is a measure of the 
direct combustion noise to the indirect combustion noise. 
Thus: 


30 


35 


40 


H,(f)=aH(f) 


( 21 ) 


H,(f)=H(f) (22) 

The direct combustion noise and the entropy noise have the 
same origin in the combustion process. Consequently it is 
assumed that the input direct combustion noise spectral den- 
sity and the input entropy fluctuation spectral density have the 


same form. 


(23) 


by: 


45 


50 


Consequently, the measured cross-spectral density is given 




(24) 


= //(/)(£->^f™ + a)C,j = 


For the frequency range of 50-250 Hz, a is negligible and 60 
H(f)=K where log^g (K) — 10. Consequently, the transfer 
function of the system is 

H,(f)=Ke-F"f‘'= (25) 

The combustor entropy noise, and combustor hydrody- 65 
namic noise, N^, are assumed independent of each other and 
independent of the tailpipe noise, Nj.. Thus: 


(26) 

A linear curve fit covering the frequency range from 
50-250 Hz was made to the cross spectrum phase angle for the 
measurement made using the combustor sensor 102 and both 
the first and second turbine exit sensor 104, 106. Specifically, 
Table 2 illustrates the linear curve fit results based on the 
results from the combustor sensor 102 and the first turbine 
exit sensor 104. Table 3 illustrates the linear curve fit results 
based on the results from the combustor sensor 102 and the 
second turbine exit sensor 106. 


TABLE 2 


Linear Fit Values 

Percent 




Std. 

Maximum 

a 

b 

Corre- 

Dev. 

Power 

ms 

Degrees 

lation 

ins 

48 

3.99 

-78.06782 

0.985 

0.0701 

54 

3.86 

-79.9248 

0.982 

0.0746 

60 

4.0229 

-90.7813 

0.974 

0.0942 

71 

3.7079 

-75.56375 

0.97 

0.093 

87 

3.4814 

-83.738 

0.948 

0.117 

98 

3.99 

-78.0678 

0.985 

0.07 

99 

3.3435 

-83.024 

0.944 

0.117 


TABLE 3 

Linear Fit Values 

Percent 

Maximum 

Power 

a 

ms 

b 

Degrees 

Corre- 

lation 

Std. 

Dev. 

ins 

48 

4.024 

-77.3482 

0.986 

0.0687 

54 

4.079 

-84.048 

0.985 

0.071 

60 

4.0664 

-83.6347 

0.986 

0.0683 

71 

4.0188 

-90.96 

0.966 

0.108 

87 

3.5794 

-83.3152 

0.961 

0.104 

98 

3.5477 

-88.01 

0.958 

0.107 

99 

3.57 

270.46 

0.966 

0.096 


The post -combustion residence time, t, expressed in ms, is 
determined from the slope of the linear curve fit expressed in 
Equation (27) below: 

P=r(360/1000)7-e (27) 


or: 

r=(P+e)/(360/1000)/' (28) 

where P is the phase angle in degrees, Q is the intercept in 
degrees, and f is the frequency. 

The post-combustion residence times are graphically 
shown in FIG. 17 for each power setting and microphone pair. 
The graph represented by the circles is for the combustion 
sensor 102 and the first turbine exit sensor 104. Similarly, the 
graph represented by the triangles is for the combustion sen- 
sor 102 and the second turbine exit sensor 106. The post- 
combustion residence time was measured as a function of 
engine power. For the example turbofan engine shown in FIG. 
1, the post-combustion residence time was approximately 4 
ms at idle and 3.4 ms at a maximum power setting. As a result, 
the measurement of post-combustion residence time has 
implications for fuel usage and system fault detection. 

The core noise components of the dual-spool turbofan 
engine were separated using coherence functions. A source 
location technique was used that adjusted the time delay 
between the combustor pressure sensor signal and the far- 
field microphone signal to maximize the coherence and 
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remove as much variation of the phase angle with frequency 
as possible. Fora 130° far-field microphone, a 90.03 ms time 
shift worked best for the frequency band from 0-200 Hz, 
while an 86.98 ms time shift worked best for the frequency 
band from 200-400 Hz. Hence, the 0-200 Hz band signal took 5 
more time than the 200-400 Hz band signal to travel the same 
distance. This suggests the 0-200 Hz coherent cross spectral 
density band is partly due to indirect combustion noise attrib- 
uted to entropy fluctuations, which travel at a low flow veloc- 
ity in the combustor until interactions with the turbine pres- to 
sure gradient produce indirect combustion noise. The signal 
in the 200-400 Hz frequency band is attributed mostly to 
direct combustion noise. The method disclosed herein is suc- 
cessful because acoustic and temperature fluctuations are 
related by a linear transfer function that includes a convective 15 
time delay. This experiment involved the measurements of 
pressure and temperature disturbances in a long tube con- 
nected to a combustor. This linear connection of entropy and 
pressure fluctuations implies the direct and indirect combus- 
tion noises are correlated at the source. 20 

Referring to FIGS. 18 and 19, a system 1800 and method 
1900 of measuring the post-combustion residence time in a 
gas-turbine engine is described respectively in accordance 
with an aspect of the innovation. The system 1800 includes a 
measurement component 1810 comprised of the combustor 25 
sensor 102 and the first and second turbine exit sensors 104, 
106, a receiving/calculation component 1820, a tabulation 
component 1830, and a computation component 1840. The 
system 1800 processes information from the sensors 
described above to determine the post-combustion residence 30 
time in the gas-turbine engine, as will be subsequently 
described. 

Specifically, at 1902, the combustor sensor 102 is placed at 
the combustor entrance, and the first and second turbine exit 
sensors 104, 106 are placed at the turbine exit, as described 35 
above in reference to FIG. 1. At 1904, the combustor sensor 
102 measures the combustor pressure at the combustor 
entrance and the first and second turbine exit sensors 1 04. 1 06 
measure turbine exit pressures at the turbine exit. At 1906, the 
cross-spectrum functions, specifically phase angle and coher- 40 
ence functions, are computed from the signals generated by 
the pressure sensors. Specifically, the receiving/calculation 
component 1810 receives the signals sent from the combustor 
sensor 102, and the first and second turbine exit sensors 104, 
106. Upon receipt of the signals, the receiving/calculation 45 
component 1810 calculates the cross-spectrum and coher- 
ence functions. At 1908, a linear curve fit is made to the 
cross-spectrum phase angle over an appropriate frequency 
range as determined by the coherence function being greater 
than 0.1. Specifically, the tabulation component 1820 deter- 50 
mines the linear curve fit and tabulates the linear curve fit 
results in a table, as illustrated in Tables 2 and 3 above. At 
1910, the tabulation component further calculates a slope of 
the linear curve fit. At 1912, the computation component 1830 
computes the post-combustion residence time, t, from the 55 
slope of the linear curve fit using the formula in Equation (28) 
above. 

Referring now to FIG. 20, a block diagram of a computer 
operable to execute the disclosed architecture is illustrated in 
accordance with an aspect of the innovation. In order to 60 
provide additional context for various aspects of the subject 
innovation, FIG. 20 and the following discussion are intended 
to provide a brief, general description of a suitable computing 
environment 2000 in which the various aspects of the inno- 
vation can be implemented. While the innovation has been 65 
described above in the general context of computer-execut- 
able instructions that may run on one or more computers. 
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those skilled in the art will recognize that the innovation also 
can be implemented in combination with other program mod- 
ules and/or as a combination of hardware and software. 

Generally, program modules include routines, programs, 
components, data structures, etc., that perform particular 
tasks or implement particular abstract data types. Moreover, 
those skilled in the art will appreciate that the inventive meth- 
ods can be practiced with other computer system configura- 
tions, including single-processor or multiprocessor computer 
systems, minicomputers, mainframe computers, as well as 
personal computers, hand-held computing devices, micro- 
processor-based or programmable consumer electronics, and 
the like, each of which can be operatively coupled to one or 
more associated devices. 

The illustrated aspects of the innovation may also be prac- 
ticed in distributed computing environments where certain 
tasks are performed by remote processing devices that are 
linked through a cormnunications network. In a distributed 
computing environment, program modules can be located in 
both local and remote memory storage devices. 

A computer typically includes a variety of computer-read- 
able media. Computer-readable media can be any available 
media that can be accessed by the computer and includes both 
volatile and nonvolatile media, removable and non-remov- 
able media. By way of example, and not limitation, computer- 
readable media can comprise computer storage media and 
communication media. Computer storage media includes 
both volatile and nonvolatile, removable and non-removable 
media implemented in any method or technology for storage 
of information such as computer-readable instructions, data 
strucmres, program modules or other data. Computer storage 
media includes, but is not limited to, RAM. ROM, EEPROM, 
flash memory or other memory technology, CD-ROM, digital 
versatile disk (DVD) or other optical disk storage, magnetic 
cassettes, magnetic tape, magnetic disk storage or other mag- 
netic storage devices, or any other medium which can be used 
to store the desired information and which can be accessed by 
the computer. 

Communication media typically embodies computer-read- 
able instnictions, data structures, program modules or other 
data in a modulated data signal such as a carrier wave or other 
transport mechanism, and includes any information delivery 
media. The term “modulated data signal” means a signal that 
has one or more of its characteristics set or changed in such a 
manner as to encode information in the signal. By way of 
example, and not limitation, communication media includes 
wired media such as a wired network or direct -wired connec- 
tion, and wireless media such as acoustic, RF, infrared and 
other wireless media. Combinations of the any of the above 
should also be included within the scope of computer-read- 
able media. 

With reference again to FIG. 20, the exemplary environ- 
ment 2000 for implementing various aspects of the innova- 
tion includes a computer 2002, the computer 2002 including 
a processing unit 2004, a system memory 2006 and a system 
bus 2008. The system bus 2008 couples system components 
including, but not limited to, the system memory 2006 to the 
processing unit 2004. The processing unit 2004 can be any of 
various commercially available processors. Dual micropro- 
cessors and other multi-processor architectures may also be 
employed as the processing unit 2004. 

The system bus 2008 can be any of several types of bus 
structure that may further interconnect to a memory bus (with 
or without a memory controller), a peripheral bus, and a local 
bus using any of a variety of commercially available bus 
architectures. The system memory 2006 includes read-only 
memory (ROM) 2010 and random access memory (RAM) 
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2012. A basic input/output system (BIOS) is stored in a non- 
volatile memory 2010 such as ROM, EPROM, EEPROM, 
which BIOS contains the basic routines that help to transfer 
information between elements within the computer 2002, 
such as during start-up. The RAM 2012 can also include a 5 
high-speed RAM such as static RAM for caching data. 

The computer 2002 further includes an internal hard disk 
drive (HDD) 2014 (e.g., EIDE, SATA), which internal hard 
disk drive 2014 may also be configured for external use in a 
suitable chassis (not shown), a magnetic floppy disk drive to 
(FDD) 2016, (e.g., to read from or write to a removable 
diskette 2018) and an optical disk drive 2020, (e.g., reading a 
CD-ROM disk 2022 or, to read from or write to other high 
capacity optical media such as the DVD). The hard disk drive 
2014, magnetic disk drive 2016 and optical disk drive 2020 15 
can be connected to the system bus 2008 by a hard disk drive 
interface 2024, a magnetic disk drive interface 2026 and an 
optical drive interface 2028, respectively. The interface 2024 
for external drive implementations includes at least one or 
both of Universal Serial Bus (USB) and IEEE 1394 interface 20 
technologies. Other external drive connection technologies 
are within contemplation of the subject innovation. 

The drives and their associated computer-readable media 
provide nonvolatile storage of data, data structures, com- 
puter-executable instructions, and so forth. For the computer 25 
2002, the drives and media accommodate the storage of any 
data in a suitable digital format. Although the description of 
computer-readable media above refers to a HDD, a remov- 
able magnetic diskette, and a removable optical media such as 
a CD or DVD, it should be appreciated by those skilled in the 30 
art that other types of media which are readable by a com- 
puter, such as zip drives, magnetic cassettes, flash memory 
cards, cartridges, and the like, may also be used in the exem- 
plary operating environment, and further, that any such media 
may contain computer-executable instructions for perform- 35 
ing the methods of the innovation. 

A number of program modules can be stored in the drives 
and RAM 2012, including an operating system 2030, one or 
more application programs 2032, other program modules 
2034 and program data 2036. All or portions of the operating 40 
system, applications, modules, and/or data can also be cached 
in the RAM 2012. It is appreciated that the innovation can be 
implemented with various commercially available operating 
systems or combinations of operating systems. 

A user can enter commands and information into the com- 45 
puter 2002 through one or more wired/wireless input devices, 
e.g., a keyboard 2038 and a pointing device, such as a mouse 
2040. Other input devices (not shown) may include a micro- 
phone, an IR remote control, a joystick, a game pad, a stylus 
pen, touch screen, or the like. These and other input devices 
are often connected to the processing unit 2004 through an 
input device interface 2042 that is coupled to the system bus 
2008, but can be connected by other interfaces, such as a 
parallel port, an IEEE 1394 serial port, a game port, a USB 
port, an IR interface, etc. 

A monitor 2044 or other type of display device is also 
connected to the system bus 2008 via an interface, such as a 
video adapter 2046. In addition to the monitor 2044, a com- 
puter typically includes other peripheral output devices (not 
shown), such as speakers, printers, etc. 

The computer 2002 may operate in a networked environ- 
ment using logical connections via wired and/or wireless 
communications to one or more remote computers, such as a 
remote computer(s) 2048. The remote computer(s) 2048 can 
be a workstation, a server computer, a router, a personal 
computer, portable computer, microprocessor-based enter- 
tainment appliance, a peer device or other common network 
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node, and typically includes many or all of the elements 
described relative to the computer 2002, although, for, pur- 
poses of brevity, only a memory/ storage device 2050 is illus- 
trated. The logical connections depicted include wired/wire- 
less connectivity to a local area network (LAN) 2052 and/or 
larger networks, e.g., a wide area network (WAN) 2054. Such 
LAN and WAN networking environments are commonplace 
in offices and companies, and facilitate enterprise-wide com- 
puter networks, such as intranets, all of which may connect to 
a global communications network, e.g., the Internet. 

When used in a LAN networking environment, the com- 
puter 2002 is connected to the local network 2052 through a 
wired and/or wireless communication network interface or 
adapter 2056. The adapter 2056 may facilitate wired or wire- 
less cormiiunication to the LAN 2052, which may also 
include a wireless access point disposed thereon for commu- 
nicating with the wireless adapter 2056. 

When used in a WAN networking environment, the com- 
puter 2002 can include a modem 2058, or is connected to a 
cormnunications server on the WAN 2054, or has other means 
for establishing cormnunications over the WAN 2054, such as 
by way of the Internet. The modem 2058, which can be 
internal or external and a wired or wireless device, is con- 
nected to the system bus 2008 via the serial port interface 
2042. In a networked environment, program modules 
depicted relative to the computer 2002, or portions thereof, 
can be stored in the remote memory/storage device 2050. It 
will be appreciated that the network connections shown are 
exemplary and other means of establishing a communications 
link between the computers can be used. 

The computer 2002 is operable to communicate with any 
wireless devices or entities operatively disposed in wireless 
communication, e.g., a printer, scanner, desktop and/or por- 
table computer, portable data assistant, communications sat- 
ellite, any piece of equipment or location associated with a 
wirelessly detectable tag (e.g., a kiosk, news stand, restroom), 
and telephone. This includes at least Wi-Fi and Bluetooth™ 
wireless technologies. Thus, the communication can be a 
predefined structure as with a conventional network or simply 
an ad hoc communication between at least two devices. 

Wi-Fi, or Wireless Fidelity, allows connection to the Inter- 
net from a couch at home, a bed in a hotel room, or a confer- 
ence room at work, without wires. Wi-Fi is a wireless tech- 
nology similar to that used in a cell phone that enables such 
devices, e.g., computers, to send and receive data indoors and 
out; anywhere within the range of a base station. Wi-Fi net- 
works use radio technologies called IEEE 802. 1 1 (a, b, g, etc.) 
to provide secure, reliable, fast wireless connectivity. A Wi-Fi 
network can be used to connect computers to each other, to 
50 the Internet, and to wired networks (which use IEEE 802 .3 or 
Ethernet). Wi-Fi networks operate in the unlicensed 2.4 and 5 
GHz radio bands, at an 11 Mbps (802.11a) or 54 Mbps 
(802.1 lb)datarate, forexample, orwithproducts that contain 
both bands (dual band), so the networks can provide real- 
55 world performance similar to the basic lOBaseT wired Eth- 
ernet networks used in many offices. 

To summarize the formation of NO„^ in gas-turbine com- 
bustors, NO,^ is formed in a distributed zone manner and that 
higher temperature-rise combustors will be required as tur- 
60 bine materials improve. While most correlation of NO„. pro- 
duction equations apply in the primary zone, research has 
found equations for other zones. The innovation disclosed 
herein is based on a combustor/turbine system study and not 
on each component separately, which does not lead to these 
65 findings. 

Others have indicated that using a constant thrust power 
setting a decrease in NO,^ emissions as a function of engine 
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age is observable. This is attributed to an increase in mass 
flow due to hot section damage. It is suggested that turbine 
damage results in lower NO^ emissions rate due to decreased 
residence time. Consequently, the post-combustion residence 
time measurement procedure disclosed herein may be uti- 5 
lized as a tool to detect turbine damage from aging or verify 
proper operation. 

The typical fault diagnostic turbine system sensor system 
depends on measuring such items as fan exit pressure, LPC 
exit pressure, burner pressure, LPC exit temperature, HPC to 
exit temperature, exhaust gas temperature, fuel flow, low 
spool speed, and high spool speed. None of these items con- 
vey the same information as the combustor residence time, 
which is a function of the turbine blade system operating 
condition and geometry. The available time to take corrective 15 
or compensatory actions such as repair or replace a part or 
reduce system operational loads to extend the life of the 
faulted part may be reduced with the additional information 
obtained from the irmovation. 

A gas turbine engine in a military or commercial aero- 20 
engine, or in industrial environment is a safety-critical sys- 
tem, which needs real-time fault detection and a decision 
support system to advise corrective actions so that the system 
can continue to function without jeopardizing the safety of 
personnel or damage to the equipment involved. Information 25 
on the status of the post-combustion residence time can pro- 
vide additional information not available from any current 
sensor used in current fault detection systems. 

Finally, in addition, jet fuel costs are 30 percent of an 
airlines cost. The status of post-combustion resident time as a 30 
function of time might provide information related to fuel 
usage. 

'\^at has been described above includes examples of the 
irmovation. It is, of course, not possible to describe every 
conceivable combination of components or methodologies 35 
for purposes of describing the subject innovation, but one of 
ordinary skill in the art may recognize that many fiirther 
combinations and permutations of the innovation are pos- 
sible. Accordingly, the innovation is intended to embrace all 
such alterations, modifications and variations that fall within 40 
the spirit and scope of the appended claims. Furthermore, to 
the extent that the term “includes” is used in either the 
detailed description or the claims, such term is intended to be 
inclusive in a manner similar to the term “comprising” as 
“comprising” is interpreted when employed as a transitional 45 
word in a claim. 

What is claimed is: 

1. A method of measuring a residence time in a gas-turbine 
engine: 

placing a plurality of pressure sensors at a combustor 50 
entrance and at a turbine exit of the gas-turbine engine; 

measuring a combustor pressure at the combustor entrance 
and a turbine exit pressure at the turbine exit; 

computing a plurality of cross-spectrum functions between 
a combustor pressure sensor signal from the measured 55 
combustor pressure and a turbine exit pressure sensor 
signal from the measured turbine exit pressure; 

applying a linear curve fit to the plurality of cross-spectrum 
functions; and 

computing a post-combustion residence time from the lin- 60 
ear curve fit. 

2. The method of claim 1, wherein prior to computing a 

post-combustion residence time from the linear curve fit, the 
method further comprising calculating a slope of the linear 
curve fit. 65 

3. The method of claim 2, wherein the slope of the linear 
curve fit is a slope of the combustor pressure sensor signal and 
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the turbine exit pressure sensor signal cross spectrum phase 
angle over a predetermined frequency range. 

4. The method of claim 3, wherein the predetermined fre- 
quency range is 0-400 Hz. 

5. The method of claim 4, wherein the combustor pressure 
sensor signal and the turbine exit pressure sensor signal is 
generated by indirect combustion noise. 

6. The method of claim 1, wherein the plurality of cross- 
spectrum functions include a phase angle function and a 
coherence function. 

7. The method of claim 1, wherein the plurality of pressure 
sensors includes a combustor pressure sensor disposed at a 
combustor entrance, a first turbine exit pressure sensor dis- 
posed at a first turbine exit, and a second turbine exit pressure 
sensor disposed at a second turbine exit. 

8. The method of claim 7 further comprising calculating a 
slope of a cross spectrum phase angle over a predetermined 
frequency range from a combustor sensor pressure signal of 
the combustor pressure sensor and a first turbine exit sensor 
pressure signal of the first turbine exit pressure sensor or a 
second mrbine exit sensor pressure signal of the second tur- 
bine exit pressure sensor. 

9. A system to measure a post-combustion residence time 
in a gas-turbine engine comprising: 

a measurement component that measures a plurality of 
combustion and/or turbine flow pressures in the engine; 

a receiving/calculation component that calculates a plural- 
ity of cross-spectrum functions based on the measured 
plurality of combustion and/or turbine flow pressures; 

a tabulation component that tabulates a linear curve fit 
based on the calculated plurality of the cross-spectrum 
functions; and 

a computation component that computes the post-combus- 
tion residence time in the gas-turbine engine. 

10. The system of claim 9, wherein the plurality of cross- 
spectrum functions include a phase angle function and a 
coherence function. 

11. The system of claim 9. wherein the tabulation compo- 
nent further calculates a slope of the linear curve fit. 

12. The system of claim 11, wherein the slope of the linear 
curve fit is a slope of a combustor pressure sensor signal from 
the measurement component and a turbine exit pressure sen- 
sor signal from the measurement component cross spectrum 
phase angle over a predetermined frequency range of 0-400 
Hz. 

13. The system of claim 12, wherein the combustor pres- 
sure sensor signal and the turbine exit pressure sensor signal 
is generated by indirect combustion noise. 

14. The system of claim 9, wherein the measurement com- 
ponent includes a combustor pressure sensor disposed at a 
combustor entrance, a first turbine exit pressure sensor dis- 
posed at a first turbine exit, and a second turbine exit pressure 
sensor disposed at a second turbine exit. 

15. The system of claim 14, wherein the tabulation com- 
ponent calculates a slope of a combustor sensor pressure 
signal from the combustor pressure sensor and a first turbine 
exit sensor pressure signal from the first turbine exit pressure 
sensor or a second turbine exit sensorpressure signal from the 
second turbine exit pressure sensor cross spectrum phase 
angle over a predetermined frequency range. 

16. A method of measuring a residence time in a gas- 
turbine engine comprising: 

placing a plurality of pressure sensors at a combustor 
entrance and at a turbine exit of the gas-turbine engine; 

measuring an indirect combustion noise of the gas-turbine 
engine; 
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computing a phase angle function and a coherence function 
of the indirect combustion noise; 

applying a linear curve fit to the plurality of cross-spectrum 
functions; and 

computing a post-combustion residence time from the lin- 5 
ear curve fit. 

17. The method of claim 16, wherein prior to computing a 

post-combustion residence time from the linear curve fit, the 
method further comprising calculating a slope of the linear 
curve fit. to 

18. The method of claim 17, wherein the slope of the linear 
curve fit is a slope of the combustor pressure sensor signal and 
the turbine exit pressure sensor signal cross spectrum phase 
angle over a predetermined frequency range. 

19. The method of claim 16, wherein measuring the indi- 15 
reel combustion noise includes measuring a combustor pres- 
sure at the combustor entrance and a turbine exit pressure at 
the turbine exit. 

20. The method of claim 16, wherein the plurality of pres- 
sure sensors includes a combustor pressure sensor disposed at 20 
a combustor entrance, a first turbine exit pressure sensor 
disposed at a first turbine exit, and a second mrbine exit 
pressure sensor disposed at a second mrbine exit. 



